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John  Pople’s  Model  Chemistries 
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The  CBS-Q//B3  Model 


Geometry:  B3LYP/6-3 1  lG(2d,d,p) 

ZPE  &  Thermal:  B3LYP/6-31  lG(2d,d,p) 

SCF:  UHF/6-3 1 1+G(3d2f,2df,2p) 

CBS(extrap)  -  MP2:  UHF/6-3 1 1+G(3d2f,2df,2p) 
MP3,  MP4(SDQ):  MP4(SDQ)/6-31+G(d(f),d,p) 
MP5  -  CCSD(T):  CCSD(T)/6-3 1+G(d’) 


Empirical  Corrections: 

-5.79  mEh  |S|2iipiC2Jii  (2e  -) 

-9.54  mEh  A(S2)  (Spin  Contam.) 

E(core)  -3.92  mEh  QNa  +  2.83  mEh  Q2Na 
Experimental  Atomic  Spin-Orbit  Interaction 
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Test  Case 

The  X  1Eg+  ground 
state,  along  with  the  a 
3IU+,  A  ’  5Zg+,  and 
A  ”  7Eu+  excited  states  of 
N2  dissociates  to  two 
4 ^3/2  ground  state  N 
atoms.  Other  low-lying 
singlet  and  triplet 
states 

dissociate  to  N  atoms 
in  2D  and  2P  excited 
states. 


R 


NN 
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Benchmarks 

The  FCI/cc-pVDZ 

energies  of  low-lying 

states  of  N2 

[H.  Larsen,  J.  Olsen, 

P  Jorgensen,  and 
O.  Christiansen, 

J.  Chem.  Phys.  113, 
6677  (2000)]  provide 
benchmarks  for 
calibration  of  model 
chemistry  candidates. 
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10.00 


The  error  in  the 
energy  of  the 
ground  state  of 
N2,  relative  to  the 
Full  Cl  energy  is 
insensitive 
to  bond  length, 
but  a  constant 
size-consistency 
correction  must 
be  made  for 
comparisons 
with  atoms. 


0) 

0) 

1— 


£ 


o 

LU 


E? 

5> 

c 

LU 

d> 

> 

oj 

d> 


N  X1!  + 
2  g 


"Approved  for  public  release;  distribution  unlimited" 


8 


Table  I.  The  effect  of  state  averaging  on  CASCF  and  CASSCF-CISD 


State 


Basis 

Set 


AE 

(mEh) 


tflg  1.10  2ZaP 

3ZaP 

A3 1*  1-27  2ZaP 

3ZaP 


CAS(6e-,6orb) 

-109.1045238  -109.0781418  26.3820 
-109.1251169  -109.0979344  27.1825 
-108.8740932  -108.8606554  13.4378 
-108.8898435  -108.8761258  13.7177 


CAS(6e.,6orb)  -  CISD 

y'Zg  1.10  2ZaP  -109.2896344  -109.2885919  1.0425 

A3Zu+  1.27  2ZaP  -109.0736612  -109.0733565  0.3047 
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Table  II.  Calculated  N2  cc-pVDZ  excitaton  energies  (mEh). 


RNN(A)  Excitation  CAS(6,6)  Qig^  FCI 


lv-  + 


1.1113  a:1! 


>  a  ‘ng 
«’  %- 


Au 

1.2700  X  %+  ^  a  ln„ 

-►  «’  %' 
— >•  w  'Au 

1.4288  X'H*  ->  a  {U„ 

-»• «’  %' 

RMS  Error 


404.8 

345.6 

343.2 

375.9 

367.9 

367.0 

412.1 

383.5 

381.7 

313.8 

274.6 

272.2 

280.6 

271.6 

270.6 

336.9 

290.3 

288.2 

252.8 

226.7 

223.4 

215.9 

207.0 

206.2 

35.3 

2.0 
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The  SCF 

energy 

converges 

to  the  CBS 

(complete 
basis  set) 
limit  as  the 
exponential 
of  the  square- 
root  of  the 
number  of 
Gaussian 
basis  functions. 
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EsCF  Limit  =  E(ll2)  + 

{exp[a(nj)1/2  -  a(n2)12]  -  l}'1  {E(n2)  -  E(itj)} 


Table  III.  Linear  SCF  extrapolation  parameters  for  the  neon  atom. 


Hi 

n2 

{ exp[a(tij)1/2  -  a(ti2)1/2]  -  l}'1 

5 

6 

0.426 

6 

7 

0.490 

7 

8 

0.551 

8 

9 

0.609 

6 

8 

0.132 

8 

10 

0.178 

10 

12 

0.220 

12 

14 

0.262 
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Table  IV.  The  convergence  of  the  N2  SCF  anisotropy  energy  (hartree) 


with  the  number  of  polarization  functions  in  the  basis  set. 


Basis  Set 

Escf 

AEscf(h) 

Escf  Error 

Extrapolation 

Factor3 

16s10p 

-108.910331 

n  HQ07Q/1 

16s10p1d 

-108.981649 

-0.071318 

-0.011416 

16s10p2d1f 

-108.991667 

-0.010018 

-0.001398 

0.139 

16s10p3d2f1g 

-108.992847 

-0.001180 

-0.000218 

0.184 

16s10p4d3f2g1h 

-108.993025 

-0.000178 

-0.000040 

0.22 

1 6s1 0p5d4f3g2h1  i 

-108.993057 

-0.000032 

-0.000008 

0.25 

SCF  Limit  = 

-108.993065 

a.  Factor  =  EScf  Error(n)  /  AEScF(n) 
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The  nZaP  BasisSets 
The  number  of 
Gaussian  primitives  in 
the  2ZaP  basis  sets  was 
selected  to  give  a 
constant,  1  millihartree, 
SCF  error  per  electron. 
The  3ZaP  basis  sets 
each  contain  two  more 
primitives  of  each 
angular  momentum. 
The  4ZaP  basis  sets 
include  an  additional 
two,  etc. 


2ZaP  Non  Relativistic  Basis  Sets 


Atomic  Number,  Z 
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The  nZaP  BasisSets 

The  contracted  nZaP  basis 
sets  are  similar  in  structure  to 
the  cc-pVnZ  basis  sets,  but 
with  diffuse  valence  functions 
added.  The  number  of 
primitives  in  each  contracted 
function  is  sufficient  to 
prevent  contraction  from 
significantly  raising  the  SCF 
energy.  All  “correlation 
primitives”  ( valence  and 
polarization)  were  optimized  at 
the  MP2  level. 


Atomic  Number,  Z 
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UHF  energies 
mimic  the 
basis  set 
convergence 

of  CASSCF 

energies,  but 

ROHF 

energies 
do  not. 


N  ^  at  R  =  3.0  A 

2  g  NN 


1  x  106  exp[-5(n  )1/2] 
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UHF’  PUHF’  CAS(10,8)’ 


Table  V.  Basis  set  convergence  of  CASSCF  and  UHFSCF  as  a  function  of  geometry. 


Rnn 

(A) 

Basis  CAS(6,6) 
Set  (Eh) 

AE 

(mEh) 

A(AE)r  uhf 

(mEh)  (Eh) 

AE 

(mEh) 

A(AE)r 

(mEh) 

1.09 

2ZaP  -109.1035124 

-108.9687667 

3ZaP  -109.1246270 

-21.1146 

0.5215  -108.9903071 

-21.5404 

0.4654 

4ZaP  -109.1279935 

-3.3665 

0.1042  -108.9935983 

-3.2912 

0.0966 

5ZaP  -109.1288265 

-0.8330 

0.0213  -108.9944093 

-0.8110 

0.0198 

6ZaP 

-108.9945297 

-0.1204 

0.0013 

1.10 

2ZaP  -109.1045238 

-108.9674580 

3ZaP  -109.1251169 

-20.5931 

0.4819  -108.9885330 

-21.0750 

0.4264 

4ZaP  -109.1283792 

-3.2623 

0.0978  -108.9917276 

-3.1946 

0.0898 

5ZaP  -109.1291909 

-0.8117 

0.0204  -108.9925188 

-0.7912 

0.0189 

6ZaP 

-108.9926379 

-0.1191 

0.0011 

1.11 

2ZaP  -109.1049575 

-108.9655455 

3ZaP  -109.1250687 

-20.1112 

0.4447  -108.9861941 

-20.6486 

0.3898 

4ZaP  -109.1282332 

-3.1645 

0.0917  -108.9892989 

-3.1048 

0.0833 

5ZaP  -109.1290245 

-0.7913 

0.0195  -108.9900712 

-0.7723 

0.0181 

6ZaP 

-108.9901892 

-0.1180 

0.0011 
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Ecas(CBS)  =  ECAS(3ZaP) 


+ 1 .2  x  [euhf  (4 ZaP)  -  Euhf  (3 ZaP)\ 


E 'cas  QZaP)  -  Ecas  (2 ZaP) 

Euhf  QZaP)  -  EUHF  (2 ZaP) 
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Table  VI.  Calculated  Full  Valence  CASSCF  energies  (hartree  au). 


2ZaP 

3ZaP 

4ZaP 

5ZaP 

6ZaP 

Limit 

c2xV 

-75.6233832 

-75.6412898 

-75.6437447 

-75.6442121 

-75.6442908 

-75.6443188 

C2  a  nu 

-75.6073350 

-75.6230710 

-75.6251470 

-75.6255495 

-75.6256266 

-75.6256541 

C2  b3Zg' 

-75.5777931 

-75.5941905 

-75.5962600 

-75.5966964 

-75.5967840 

-75.5968152 

CjA'n. 

-75.5627890 

-75.5788475 

-75.5809697 

-75.5813798 

-75.5814615 

-75.5814906 

CN  X2S+ 

-92.3544025 

-92.3730997 

-92.3761806 

-92.3766716 

-92.3767524 

-92.3767812 

co  xV 

-112.8946223 

-112.9181416 

-112.9218328 

-112.9224855 

-112.9226137 

-112.9226594 

n2  x%+ 

-109.1166132 

-109.1373292 

-109.1406653 

-109.1414826 

-109.1415931 

-109.1416324 

NO  X2IIi 

-129.3892421 

-129.4106758 

-129.4138042 

-129.4143983 

-129.4144941 

-129.4145282 

o2  x3zg 

-149.7411441 

-149.7681917 

-149.7713131 

-149.7719475 

-149.7720640 

-149.7721055 

02  aAg 

-149.7059304 

-149.7327122 

-149.7358092 

-149.7364403 

-149.7365522 

-149.7365920 

03  X'A, 

-224.5478659 

-224.5938866 

-224.5991383 

-224.6002546 

-224.6004950 

-224.6005806 

oh  x2n 

-75.4325178 

-75.4440618 

-75.4456305 

-75.4458720 

-75.4459251 

-75.4459441 

OH  A2E+ 

-75.2670507 

-75.2808979 

-75.2825953 

-75.2828559 

-75.2829216 

-75.2829450 

RMS  Error  0.0224147 

0.0030328 

0.0005829 

0.0001180 

0.0000310 
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Table  VII.  Extrapolated  Full  Valence  CASSCF  energies  (hartree  au). 


3ZaP 

Error 

4ZaP 

5ZaP 

C2X1Ig+(Re) 

-0.0001705 

0.0000255 

-0.0000059 

C2  a  nu(Re) 

-0.0001753 

0.0000582 

0.0000012 

C2  b3Sg(Re) 

-0.0002094 

0.0000064 

-0.0000040 

C2A1nu(Re) 

-0.0000459 

0.0000176 

0.0000026 

CN  X2E+(Re) 

-0.0001536 

0.0000366 

-0.0000031 

N2  X1Sg+(Re) 

-0.0005054 

-0.0001428 

-0.0000276 

NO  X2n,(Re) 

-0.0001124 

-0.0001120 

-0.0000175 

02  X3Ig'(Re) 

0.0000564 

-0.0001285 

-0.0000186 

02  a^gfRe) 

-0.0000230 

-0.0001280 

-0.0000168 

OH  X2n(Re) 

-0.0000524 

-0.0000353 

-0.0000136 

OH  A2S+(Re) 

0.0000168 

-0.0000603 

-0.0000194 

RMS  Error 

0.0002327 

0.0000956 

0.0000245 

Before  Extrap 

0.0030328 

0.0005829 

0.0001180 
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MP2  CBS  Limit 

The  aa  second- 
order  Moller-  0 
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MP2  CBS  Limit 

The  a(3  second- 

order  Moller- 

0) 

Plessett  electron 

0) 

1— 

t 

m 

pair  correlation 

vv 

r 

energies  converge 

to  the  CBS  limit  as 

TO 

the  inverse  third 

V1 

power  of  the 

s 

/I  \ 

maximum  angular 

<D 

Ctf5 

momentum 

01 
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Convergence  of  N2  a.p  Pair  Energies 
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Table  VIII.  The  basis  set  convergence  of  dynamic  correlaton  energy. 


State 

Basis 

Set 

e<2) 

(mEh) 

ae(2) 

(mEh) 

Ecas-cisd  - 
Ecas  (mEh) 

A(Ecas-cisd  —  Ecas) 

(mEh) 

2ZaP 

-316.0189 

-181.9089 

1.1  OA 

3ZaP 

-379.3504 

-63.3315 

-244.4998 

-62.5909 

4ZaP 

-403.7559 

5ZaP 

-413.3317 

6ZaP 

-417.4648 

CBS  Limit 

-426.90 

2ZaP 

-250.9989 

-195.6356 

1.27A 

3ZaP 

-310.1058 

-59.1069 

-255.0882 

-59.4526 

4ZaP 

-332.4087 

5ZaP 

-341.1121 

6ZaP 

-344.7652 

CBS  Limit 

-353.59 
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E'  * <■  ^(CBS)  =  ECAS_asn(3ZaP) 


CAS-CISD 


+ 


E{2)(CBS)-E{2)(3ZaP )] 


x 


'CAS-CISD 


(3ZaP)  -  E 


CAS-CISD 


(2  ZaP) 


E(2)  (3ZaP) -  E{Z>  (2ZaP) 


(2) 
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Table  IX.  Calculated  CAS-CISD  geometries  (Re  in  A)  for  low- 
lying  excited  states  of  N2. 


State 

2ZaP 

3ZaP 

CBS  limit 

Exp. 

X A+ 

1.1175 

1.1039 

1.0991 

1.0977 

A3ZU+ 

1.3066 

1.2951 

1.2905 

1 .2866 

B3ng 

1.2314 

1.2201 

1.2152 

1.2126 

w3au 

1 .2972 

1 .2859 

1.2815 

B’32u 

1 .2955 

1 .2845 

1.2800 

1 .2784 

a  > 1  Zu- 

1.2914 

1.2805 

1.2759 

1 .2755 

a  1  ng 

1 .2390 

1.2271 

1.2213 

1 .2203 

1 .2837 

1.2728 

1 .2684 

1.268 

RMS  Error 

0.0181 

0.0066 

0.0020 
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Table  X.  Calculated  CAS-CISD  harmonic  vibrational  frequencies 
(coe  in  cm"1)  for  low-lying  excited  states  of  N2. 


State 

2ZaP 

3ZaP 

CBS  limit 

Exp. 

2319.1 

2345.2 

2358.1 

2358.57 

a  3  z: 

1418.5 

1439.7 

1450.8 

1460.64 

b  3  ng 

1707.4 

1715.3 

1727.4 

1733.39 

w3a„ 

1482.0 

1499.1 

1508.9 

1501.4 

B’3ZU- 

1494.5 

1509.9 

1518.9 

1516.88 

a  >  ‘Zu- 

1508.6 

1528.4 

1538.8 

1530.25 

a  1IJg 

1648.2 

1686.1 

1696.4 

1694.21 

1543.5 

1561.3 

1571.3 

1559.26 

RMS  Error 

31.1 

12.1 

7.3 
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Table  XI.  Calculated  CAS-CISD  excitation  energies  (Ev)  of  N2. 


State 

2ZaP 

3ZaP 

CBS  limit 

Exp. 

A3ZU+ 

5.843 

6.093 

6.192 

6.224 

B3ng 

7.275 

7.397 

7.447 

7.392 

w3au 

7.233 

7.364 

7.432 

7.415 

B’3ZU- 

8.067 

8.230 

8.265 

8.217 

a  >  'lu- 

8.427 

8.471 

8.407 

8.450 

a 

8.547 

8.653 

8.562 

8.590 

8.932 

9.008 

8.927 

8.939 

RMS  Error 

0.176 

0.065 

0.037 
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Table  XII.  Calculated  CAS-CISD  dissociation  energies  (Ev)  of  N2. 


State 

2ZaP 

3ZaP 

CBS  limit 

Exp. 

A3ZU+ 

5.857 

6.093 

6.192 

6.224 

B3ng 

7.256 

7.397 

7.447 

7.392 

w3au 

7.301 

7.364 

7.432 

7.415 

B’3ZU- 

8.063 

8.230 

8.265 

8.217 

a  >  'lu- 

8.417 

8.471 

8.407 

8.450 

a 

8.537 

8.653 

8.562 

8.590 

8.914 

9.008 

8.927 

8.939 

RMS  Error 

0.167 

0.065 

0.037 
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